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Abstract. 

We investigate the motion, near the equilibrium configurations, ol an initially spinless rigid 
body subject to an external tidal field. Two cases are considered: when the center of mass of the 
body is at rest at the equilibrium point of the field generated by a generic mass distribution, and 
when it is placed on a circular orbit subject to a spherically symmetric potential. A complete 
analysis of the equilibrium configurations is carried out for both cases. First, we derive the 
conditions for the stable equilibria, and then we analyze the frequencies of oscillations around the 
equilibrium positions. In view of these results, we consider the problem of alignment of galaxies 
in clusters. After estimating the period of the oscillations induced on the galaxies by the tidal 
field of the cluster, we discuss the possible effect of resonances between stellar orbits inside the 
galaxy and the oscillations of the galaxy as a whole; this may be a mechanism responsible for 
producing an intracluster stellar population. 
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1 Introduction 

Informations about the structure and evolution of elliptical galaxies (Es) are com- 
monly inferred from their morphology. Roughly speaking, this is justified by the 
fact that their internal two-body relaxation time is longer than the Hubble time, 
so that the shape and internal dynamics of such systems probably reflect the con- 
ditions present immediately after the end of their formation process (Binney & 
Tremaine 1987, hereafter BT87). However, a problem with this scenario is that Es 
in clusters are not isolated objects, being subject to the effects of their environ- 
ment. Since these effects are acting on time scales shorter than the Hubble time, 
it is of the utmost importance to be able to isolate them when using the present 
morphology to infer the conditions at the galaxy formation epoch. In particular, 
the cluster's tidal field (TF) acts continuously, and affects in general any galaxy 
in the cluster. Thus, although the TF is not as strong as other interactions (e.g., 
galaxy-galaxy encounters), its effects can become important on the long term. 

A first possibile effect of the TF is a modification of the shape of galaxies 
belonging to the cluster itself. Ciotti & Dutta (1994, hereafter CD94), using N- 
body numerical simulations, have shown that this is not the case, since Es behave 
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like rigid bodies in the cluster TF, at variance with disk galaxies (Valluri 1993). 
Such a conclusion seems to be confirmed by observational results, because no 
systematic differences are found between cluster and field Es (for a review see 
CD94 and references therein). 

A second possible effect on galaxies is their alignment with the cluster major 
axis, or with the radial direction in the case of a nearly spherical cluster. To 
this respect, Hawley & Peebles (1975) reported "a possible indication that the 
galaxies are preferentially aligned along the radius vector to the center of the 
[Coma] cluster"; this result was confirmed by Thompson (1976). MacGillivray 
& Dodd ( 1979a, b) observed an analogous alignment for Es in the cluster centered 
on NGC439; in this case they found also a number of galaxies with their major axis 
perpendicular to the radial direction of the cluster. The same authors, observing 
another spherical cluster, found a preferential alignment of Es' major axes with the 
cluster's radial direction. Adams, Strom & Strom (1980) found a general trend for 
Es to be aligned with the cluster's major axis in 7 very elongated clusters. They 
also reported "a small but significant number of Es with their major axes perpen- 
dicular to the cluster's major axis". Fong, Stevenson & Shanks (1990) found a 
moderately significant radial alignment of Es, together with a less significant ten- 
dency of Es to be oriented perpendicularly to this direction. The best evidence for 
alignment is for the Brightest Cluster Ellipticals: Trevese, Cirimele & Flin (1992) 
described a "strong alignment of the brightest galaxy major axis with the long axis 
of the parent cluster". 

A third effect strictly related to the previous one is the possible resonant inter- 
action between stellar orbits inside the galaxies and the cluster TF: if this is the 
case then we can speculate on the possibility of collisionless evaporation of stars 
from galaxies in clusters. 

The N-body simulations of prolate Es in radial and circular orbits inside clusters 
show that an alignemnt, due to the cluster TF, can in fact be expected (CD94). 
In this paper we establish on more firm grounds the analytical discussion of the 
results of CD94 concerning the alignment. In particular, we discuss the equilibrium 
positions of a triaxial galaxy at the center of a triaxial cluster, and of a triaxial 
galaxy in circular orbit inside a spherical cluster, analyzing also the stability of 
these equilibrium configurations. For small displacements near the stable positions 
we obtain the oscillation frequencies, which are then compared with the internal 
orbital periods of the stars. Throughout our discussion (following the N-body 
results of CD94), a galaxy is described as a rigid body (RB) placed in an external 
gravitational field. By linearizing the external field, the problem is reduced to 
that of the integration of the equations of motion for a RB in a tensorial force 
field, corresponding to the TF. Although this problem is not integrable, one can 
still determine the equilibrium positions and their stability, along with the proper 
frequencies of the system. 

The paper is organized as follows. In Section 2 we assume that the center of 
mass of the RB is at rest inside an unspecified external density distribution. In 
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Section 3 we study the more complicate case of a RB on a circular orbit inside 
an external, spherically symmetric potential. In both cases, we first determine 
the equilibrium configurations, and then we study their stability. For the stable 
positions the frequencies are determined analitically. Finally, in Section 4 we 
discuss the astronomical application of the previous analysis to the problem of the 
alignment of galaxies inside galaxy clusters, and of galaxies on circular orbits within 
spherical clusters. In addition, we compare the mean stellar orbital times inside 
galaxies with the frequencies induced by the TF, and we consider the possibility 
of a collisionless evaporation induced by resonances. 



2 Body at rest at the equilibrium point of an external potential 

We start discussing the simple case of an initially spinless RB with its center 
of mass at rest in an equilibrium point of an external potential. We introduce 
an inertial reference system C, the time-independent density distribution p(x) 
producing the external field g = — V<^(x). 

Without loss of generality we can assume that the equilibrium point of the 
potential is at x = 0. We also assume that the characteristic dimensions of the RB 
are much smaller than the typical length scale over which the external potential 
changes. Thus, near the equilibrium position 

0(x) = 0(O)-i(Tx,x) + O(||x|| 3 ), (1) 

where the brackets indicates the usual inner product in K 3 and 

a 2 0( x ) 



dxidxj 



(2) 



is the TF tensor. Note the minus sign in eqs. (|l]) - (Q), due to the fact that the 
TF is defined as the truncation to the first order of the series expansion of g. We 
call prb(x;£) the density distribution of the RB of total mass M, and we assume 
that the RB center of mass Rcm is initially at rest at x = 0. The acceleration of 
the center of mass is given by 

MA CM = J gPRBd 3 x a Mg(R CM ) + T (R CM ) / (x - R C M)PR B d 3 x. (3) 

It is evident that the last integral vanishes, i.e., the motion of Rcm is not perturbed 
in the TF approximation; it is only in higher order terms of the force expansion 
that the coupling appears. As a result of the initial conditions imposed, only 
rotational motions around x = are possible. As usual we indicate with 

%j(t) = J (||x|| 2 (% - XiXj)pR B (x;i)(f 3 x, (4) 

the generic component of the RB inertia tensor Q, as seen from C. 
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2.1 Equilibrium Configurations 

In order to find the equilibrium configurations of the RB in the TF we make use 
of the second cardinal equation of Dynamics, and we search for the vanishing of 
the torque due to the TF. The tidal torque N is obtained by the series expansion 
of its integral definition: 

Ni = -e ijk J x j d Q^ /)RB(x;t)d 3 x, (5) 

where ey^ is the Levi-Civita tensor (here and in the following, the summation 
over repeated indices is assumed). Up to now no assumptions have been made 
about the orientation of the TF with respect to C. Without loss of generality 
we can assume that T is in its diagonalized form, so that, using eq. (|) and the 
antisymmetry of e^, eq. (||) can be rewritten as 

3 

Ni = Yl ^%fcAT jfc , (6) 
j,k=i 

where AT is the antisymmetric tensor of components ATjj = Tj — Tj, and Tj are 
the three principal values of the TF. 

From eq. (|6|) it is evident that the torque along a particular axis is affected only 
by the difference between the TF components along the two orthogonal directions. 
It is also obvious that if the external gravitational field is spherically symmetric 
around x = 0, or if the RB itself is spherically symmetric, then no net momentum 
acts on the RB, i.e., all its positions are of equilibrium. 

We use eq. @ to find the equilibrium configurations of the RB. Let us define 
an orthogonal reference system C with its axes directed along the principal axes 
of inertia of the RB, and introduce the orthogonal transformation matrix 0{t) 
between C and C. In C the RB is at rest, and its inertia tensor is diagonal, with 
the three principal components (Ji, 12,13). We consider only a non-degenerate 
situation (i.e., I\ 7^ I2 7^ -^3 and T\ 7^ T2 7^ T3). If O acts on the vectors £ of C as 
x = 0$,, then the following relation holds: 



% fe (i) = ]T O jfl {t)O k ,{t)I,. (7) 
n=l 

From eq. (0), assuming ATj k 7^ 0, we conclude that the RB is in equilibrium if 
and only if = for j 7^ k, i.e., if and only if 3 s is diagonal. This is true if and 
only if C and C are coincident (apart from a renaming of the RB principal axes 
of inertia). Thus, all the equilibrium positions correspond to the RB inertia axes 
oriented along the TF principal directions. 
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2.2 Stability 

We now consider the Lagrangian C = T — U associated to our problem. Expanding 
the external potential as in eq. (Q), from eq. (0) and eq. (^) one obtains 

U = J ^prbCxj^x ~ M0(O) + 7tGp(0)Trm + ^(i)/^, (8) 

where TV [3 s ] is the trace of the RB inertia tensor, and the property of diagonality of 
T in C has been exploited. The first two terms in the r.h.s. of eq. (§) are constant 
and from now on they will be dropped. Before we proceed, it is necessary to write 
explicitely the transformation matrix O between C and C . Instead of using the 
usual representation of O in terms of the Euler angles (i.e., the choice of the 3-1-3 
rotations)[] we adopt the counter-clockwise rotations 1-2-3, namely 

= O x {<p)0 2 (#)O z W), (9) 

with 

/I \ 

Oi(tp) = cos <p —sin (p , (10) 
\ sin ip cos p J 

(cos i? sin •& \ 
10, (11) 
— sin i? cos i? / 

(cos -0 — sin tp \ 
sinV cos^ . (12) 
l) 

The angular velocities associated with the rotations around the three inertia 
axes are u>J = (0,0,0), etc, where the suffix 'T' means 'transpose'. By vector 
composition, the angular velocity of the RB, as seen from C", is: 

w = ^ + Oj(V)K + Oj(^)^]. (13) 

The kinetic energy of a RB with its center of mass kept fixed is given in C by 

2T = ii(i?sin?/> + (p cos $ cos if)) 2 + I2 ("!? cos ^ — 0cosi?sin'i^) 2 + Is(ip + (psiwd) 2 . 

Similarly, the potential energy U can be expressed as: 

2U = (cos (p cos ip — sin ip sin 1? sin ip) 2 AT21 A/21 + sin p 2 cos i? 2 AT21 A/31 + 
(sin</?cos V' + cos</?sin7?sin?/j) 2 AT3iA/2i + cos ip 2 cos $ 2 AT3i A/31. 

(14) 



1 This choice complicates the treatment due to the indeterminacy of the angles ip and ip for a 
null inclination ($ = 0). 
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where we have introduced the antisymmetric tensor A y - = Jj — Ij, and neglected 
an additive constant. 

As proved in §2.1, apart from a renaming of the inertia axes, the equilibrium 
positions are given by (tp, ip) eq = (0,0,0). The linearized kinetic and potential 
energy near this point are simply 

T ~ l -h^ 2 + + \htf, (15) 

U ~ -iAr 32 A/ 3 2¥' 2 - ^AT 31 AI 31 tf 2 - iAT 21 AI 21 ^ 2 . (16) 

Note how both the quadratic forms are already diagonal: this means that the 
chosen coordinates (angles) are also the normal coordinates for the problem. The 
differential equations of the motion are readily obtained: 

.. Ar 32 A/ 3 2 n _v 

ip = (p (17) 

h 

9 = AT31 T Ahl # (18) 

V = j ^, (19) 

and their solution is trivial. Note that in the determination of stability positions 
and relative frequencies, only the ratios of the inertia moments matter. Accord- 
ingly, in order to simplify our successive discussion in §3, we define 

U =V V = T 3 - (20) 

Without loss of generality, let us assume that T% > T 2 > T 3 , i.e., that AT 32 , AT 3 i, 
and AT 2 i are all less or equal to zero. Thus, from eqs. (17)-(^), stable motions are 
possible if and only if 1 > v > u. In conclusion, a stable equilibrium configuration 
(in the non-degenerate case) is achieved when: 

1. The principal axes of the the TF are superimposed with those of the RB; 



2. If Ti > T 2 > T 3 then h < h < h- 



For what concern the behavior of the frequencies, it is obvious that their values 
increase if we increase the differences between the TF eigenvalues Ti and/or the 
values of the inertia parameters u, v. In §4 the physical meaning of points 1 and 2 
above will be discussed for our astronomical application, and an explicit calculation 
of some typical frequencies will be derived for galaxy and clusters models matching 
observational constraints. 
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Fig. 1. The three reference systems used in the discussion of the RB in circular motion inside a 
spherical potential. 

3 Body on a circular orbit in a spherical potential 

In this section we discuss the case of a triaxial RB on a circular orbit inside a 
spherically symmetric potential. Note that the assumption that the RB's orbit 
remains circular is justified by the fact that the external field, truncated at the 
TF approximation, does not affect the motion of Rcm (see §2). 

Without loss of generality, we can assume that the orbit of Rcm is in the (x, y) 
plane of the inertial system Co, the external density is given by p = p(r) and its 
associated mass inside the radius r is M(r). The natural frame for our problem 
(Fig. 1) is the non-inertial reference system C = (0;ei,e2,e3) with its origin 
placed at Rcm 5 e i = Rcm/||Rcm||) and e3 parallel to the z-axis of Cq. As a 
consequence C is rotating with the same angular velocity fl of Rcm : 

fl T = (o,o, n)> 

We indicate with £ and £ respectively the generic position and velocity vector in 
C, related to x and x by the usual transformation rules involving the orthogonal 
transformation matrix C between Co and C (see, e.g., Arnold 1978). 

3.1 Equilibrium configurations 

As in the previous case we start discussing the RB equilibrium positions, by means 
of the second cardinal equation of Dynamics. The relations between the angular 
momenta and torques in the two systems are given by 

L = MR CM A V CM + C(L + 9tt), (22) 

and 

N = MR CM A A CM + CN, (23) 



r= Rpm 



(21) 
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where 9 is the (time-dependent) RB inertia tensor in C, and 

L = J (i A kpKBdh, N = J I A (C T g)p RB d 3 H. (24) 

From Lo = No one obtains, upon multiplication by C T : 

L = N — 9 17 — 917 — 17 A L — 17 A 917. (25) 

The equilibrium positions are determined by the request that L = L = 9 = 0, 
along with the condition (verified by hypothesis) that 17 = 0. Eq. ( p5|) then gives 



N-ttA9tt = 0. (26) 

The first term to be discussed is the tidal torque N. A treatment analogous to 
that described in §2 gives again eq. (pi), where now (CD94) 



3q-2 s 

T = -Z^JL | o 1 | , (27) 
1 



A-nGp 



having defined 



where 



qi r) - g| < I 2 , (28) 

3M(r) _ 30V) 
Plrj - 4^ " 4vrG ' 1 9j 

Note how in this case the TF is degenerate, i.e., two of three TF principal 
components are equal: T2 = T3 and so AT23 = 0; also, AT12 = AT13. 
We can now proceed with the discussion of eq. (f26|), which gives 



9 23 ^ 2 = 0, (30) 

9i 3 (AT 13 + ft 2 ) = 0, (31) 

9i 2 AT 12 = 0. (32) 

Since from eqs. (|27| ) and (2) we get 

AT13 = AT 12 = 3(1 - q)Q 2 > 0, (33) 



eqs. (30) - ( p2[) give the same conditions as in the previous section, namely the 
equilibrium is achieved when the off-diagonal components of 9 vanish in C, which 
in this case means that one inertia axis is directed radially, and another is directed 
along the angular velocity direction. Thus, 17 removes the intrinsic degeneracy of 
the TF. 



2 If dp/dr < then the inequality q < 1 is trivial to prove. If p(r) oc r a (a < 3), then 
q = 1 — a/3. Outside a spherical body q — 0. 
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3.2 Stability 

We will now linearize C around the equilibrium points. Let C be the reference 
system whose axes are the principal axes of the RB, and O the transformation 
matrix between C and C' [see Fig. 1 and eq. @]. The final expression for the 
kinetic energy is: 

T = MVSu , (0 + + (u) 

2 2 ' 1 ; 



where w, given by eq. (|l3|), is the angular velocity of the RB as seen from C", and 
f2' = T f2. 

Obviously, the first term in the r.h.s. of the previous equation is constant, and so 
can be discarded. Equation ( |34|) contains a term linear in the generalized velocities, 
a well known consequence of the time dependence of the new coordinate system. 
With a careful analysis, however, it is easy to prove that this linear component 
is identically null for a constant 17, so that it does not affect the motion. The 
potential energy remains unaffected by the coordinate change, and its expression 
is again given by eq. (|0]) . The differential equations associated with the quadratic 
part of C are: 

if = — <p+ , (35) 

» = (o 2 + Ar 13 )A/3i ^ _ n(i 1 + i 2 -i 3 )<p 

ip = V; (37) 

note that for = we reobtain eqs. (|l7|)-(|l9|) (with AT23 = 0). 

We now analyze the stability of the system described by eqs. (|35|)-([37|)P|. The 
first result follows directly from eq. (|37|), from which, by virtue of the assumption 
AT 2 i < 0, we get the stability condition I 2 > I\. 

Equations (|35|)-(36), after some standard manipulations, are separated into two 



linear fourth-order equations with constant coefficients, that are solved setting 

ip(t) = Xexp(Xt), -&(t) = Yexp(Xt). (38) 
The resulting characteristic equation for A is biquadratic^: 

A 4 + a\ 2 + = 0, (39) 

3 These linearized equations can also be obtained by expressing eq. ( |25| ) in C": 

Iu = N' - {ill + us) A 7(17 +u)- I{S7 A w). 

Note the identical formal aspect of this system with that discussed in BT87, p. 138 for the stability 
of Lagrangian points. 

4 The separation introduces two additional integration constants for each variable. These are 
determined by the requirement that eqs. (j3q)-(p6[) are satisfied at each instant of time. 
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where 



a = Q 



(u + v- l) 2 + 1-v + (4-3g)(l-M) 



uv u 



(40) 



^ (4 -3g)(l -«)(!-«) 



/? = ^ — (41) 

uv 

and u and i> are given in eq. ( p0|) . Any solution of eq. (^) is of the form 

\ = n\(u,v;q), (42) 

where A is a dimensionless function, and the dependence on r is contained in Q and 
q. This means that in the determination of stable equilibria only (u, v; q) matter, 
i.e., the parameter space is three-dimensional. Moreover, for fixed (u,v;q) the 
frequencies are proportional to 0, the orbital velocity of Rcm around the center 
of the external density distribution. 

Equation ( p^ ) has four solutions: if A is a root also —A is a root, and so the 
general solution will contain an exponential increasing term if 3Re(A) ^ 0, making it 
unstable. The only possibility for a stable equilibrium is then for the four solutions 
to be purely imaginary, and this is possible if and only if the two solutions of the 
associated quadratic equation are real and negative. The reality requires that 

A = a 2 - 4(3 > 0, (43) 

and for the negativity it is required that 

a > 0; p > 0. (44) 



Note that the conditions (j43[)-(f44[) are formally equivalent to the single condition 
a > 2^/]3. The positivity of (3 is obtained in two different cases: 

(a) u < v < 1, i.e., I3 > I2 > h- 

(b) 1 < u < v, i.e., I 2 > h > h- 

where the stability of tp motion (I2 > I±) has been already assumed. 

Case (a) is of easy discussion: in fact, a is a sum of positive quantities, and A 
can be rewritten as a sum of positive quantities as well, independently of the value 
of q{r). Thus, the motion around the equilibrium position is stable for I3 > I2 > l\, 
for any given density profile of the external density distribution. This corresponds 
to region (a) in Fig. 2. 

Case (b) is far more interesting. Its discussion is algebrically cumbersome, but 
not intrinsecally difficult. It can be proved that when 1 < u < v the request A > 
implies a > 0, so that the only condition to be studied is A > 0. The stability 
region in the (u, v) plane is plotted in Fig. 2 for a specific value of q. Note that 
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Fig. 2. The (it, v) plane with the various critical curves, for q — 0.9. In region a the major axis 
of the RB is directed toward the center of the external potential, the intermediate axis is directed 
tangentially to the orbit and the minor axis is perpendicular to the orbital plane (Fig. 1). In 
region b = bi U 62 the major axis of the RB is perpendicular to the orbital plane, the medium axis 
is directed toward the center, and the short axis is parallel to the orbital velocity. 

the stability region relative to case (b) is made of two qualitatively different parts, 
separated by a vertical asymptote, located at: 



In region b± (1 < u < u a , see Fig. 2) there is no upper bound on I2/I3, and 
q — > 1 =>■ u a — > 00, i.e., region b\ expands to all v > u > 1. For any finite value of 
q, there exist another stability region (62, see Fig. 2) to the right of u a , whose lower 
bound is again represented by v = u. The upper bound on I2/I3 (determined by 
the condition A = 0) is given by 



The functions u a (q) and u c (q) are plotted in Fig. 3. A very important case is 
obtained for q = 0, i.e., in the case of a RB orbiting outside a spherical body: in 
this case the region of equilibrium is very small, since it a (0) = 1/2 + 1/^- 1.08 
and n c (0) = 1/2 + 3y/E/10 ~ 1.17. 

3.3 Frequencies 

Having found the criteria for stable equilibrium, we now focus our attention on 
the frequencies of the motions around the (ip, ip) axes. The frequency for the t/j 
motion is obtained directly from eq. ( fl9| ) as 




(45) 




(46) 




(47) 
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Fig. 3. The behaviour of M a and u c (see Fig. 2) as a function of q. 

^ = nJk(v-u); k = 3(1 - g). (48) 

Its value increases by increasing f2 and/or decreasing q, i.e., by making the external 
density more concentrated. 

The motion around ip and $ results from the superposition of two independent 
oscillations, with frequencies 

2n_ = 2vr_ = ^ / l + (<r + fc)<r ± ^/l + <t*(t + kf - 2a{r + 2kr - k) 
where 

1 — u 1 — v 
a = ; r = . 

v u 

We stress again the presence of the proportionality constant O for all frequen- 
cies: assuming that the dimensionless factor is of the order of unity, the oscillatory 
motion around the stable positions is characterized by periods of the order of the 
RB orbital time. 

4 An Astrophysical Application: Galaxies in Clusters 

We now apply the previous results to the cases of a triaxial galaxy at the centre 
of a triaxial cluster, and of a triaxial galaxy in circular orbit inside a spherical 
cluster. We start presenting the galaxy model used. 

4.1 The galaxy model 

We assume that the galaxy density is stratified on homeoidal surfaces labeled by 

3 x 2 

m 2 = ^-|, ai > « 2 > Q 3 , (50) 
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i.e., p g = p g (m). In the reference system implicitely assumed in eq. (^), all the 
non-diagonal terms of 9 vanish, and the diagonal terms are given by: 

k = —aia 2 a 3 (a 2 j + al)h g ; i / j / k, (51) 
where the form factor is 

poo 

h g = pg(m)m 4: dm. (52) 
Jo 

Note that the inertia moments satisfy I\ < I 2 < ^3- Note also that in all our 
results in §2 and §3 only the quantities u = I1/I2 and v = I2/I3 appear, and in 
particular we conclude that for a triaxial homeoidal body the equilibrium positions 
and the relative frequencies do not depend on the quantity h g , i.e., they do not 
depend on the particular density distribution under consideration. 

However, a particular density distribution for the galaxy model has to be 
assumed in order to estimate the characteristic stellar orbital times, and so we 
adopt an ellipsoidal generalization of the widely used 7-models (Denhen 1993, 
Tremaine et al. 1994): 

PgM = aia 2 a 3 ^m^l + m)^' (53) 

where < 7 < 3. The density in the outer parts goes like m~ 4 independently 
of 7; for 7 = the models present a flat density core. It is straightforward 
to show that a density stratified on constant m surfaces, when projected on its 
principal planes, originates a brightness distribution with elliptical isophotes whose 
ellipticity is the same as that of the corresponding section of the ellipsoid. For 
example, projecting along X3 the isophotal axial ratio in the x\ — x 2 plane is 
012/011, and the corresponding ellipticity is €21 = 1 — 02/01 • This property is useful 
when constructing realistic galaxy models: mean observed values are e = 0.2 — 0.4 
(corresponding to E2-E4 galaxies). 

For the models described by eq. Q53|) the mass inside m is given by: 



Mg(m) = 4 TOl a 2 a 3 / p g (t)t 2 dt = Mg (54) 
Jo (1 + mf 7 



= T^ZZITZ ^M I ^3- 7 ' ( 55 ) 



the mean density inside m is 

aia2«3 47r m 7 (l + m) 



and an estimate of the order of magnitude of the characteristic orbital time of stars 
inside m is finally obtained as usual as 



/ o7T 

Poib(m) ~ 4P dyn (m) = J -, r , (56) 
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Fig. 4. Extimated orbital times as function of m for various 7. The curves from top to bottom 
correspond, respectively, to 7 equal to 0, 1, 1.5, 2, and 2.5. The black diamonds correspond to the 
half-mass orbital time, i.e., the orbital time calculated at mi/2 = p 1 ^ 3-7 ' — 1] _1 , the homeoid 
containing half of the galaxy mass. 



(BT87, p.37)Q. 

From eqs. (|55|)-(l56|) one can write 



-forb = \ / orb — 3.0 



a »° (1 - £2l)(1 -^P orb(m , 7 ) xl oSr, (57) 



^,11 

where M gj n is the galaxy mass normalized to 10 11 M Q and 0:1,10 is the galaxy core 
major axis in 10 kpc. In Fig. 4 we plot the dimensionless function -P orD , for various 
values of 7. 

4.2 Triaxial Galaxies at the Center of Triaxial Clusters 

We now focus on the TF produced by a triaxial cluster. We assume that the 
cluster's density is stratified on homeoidal surfaces p c = p c (^ 2 ), with 

3 x 2 

^ 2 =E-|; ai >a 2 >a 3 . (58) 

The associated potential is found using the well-known formulas from potential 
theory (Kellog 1953; Chandrasekhar 1969, Cap.3; BT87 p.61). With the aid of the 
auxiliary function 

^{l 2 ) = / p c (t)dt, (59) 
Jo 

5 Note that, strictly speaking, the previous formula applies only to the case of a point mass 
in radial or circular orbit inside a homogeneous, spherically symmetric density distribution. As 
pointed out by the referee, particular attention should be paid in the case of an incompressible, 
self-gravitating, prolate spheroid steadily rotating about its axis of symmetry: in fact, in this case 
the equations of hydrodynamics do not admit solutions (see, e.g., Florides & Spyrou 1993). 
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the potential can be written as 

Jo A(t) 

where 

A(r) = y/(al + T)(a% + T)(4 + T), (61) 

and 



T.*) = ±-fc. (62) 



The fact that the center of the cluster is an equilibrium point is trivial to prove, 
and with very simple algebra the following integral expression for the principal 
components of the TF at the cluster center are found: 

f°° dr 

Ti = -2irGp c (0)aia 2 a 3 / — 5— — . . . . (63) 







(a? + r)A(r) 



Note that the particular form of the cluster density distribution does not enter 
eq. (p^l), and that if a\ > a 2 > (13 then T\ > T 2 > T3, as assumed in §2.2. We 
can now give the astrophysical interpretation of the results obtained in §2.2: the 
stable equilibrium position of a triaxial galaxy at the center of a triaxial cluster 
corresponds geometrically to the three principal axes of the galaxy to be collinear 
- in the same order - with the principal axes of the cluster density distribution. 
This configuration reproduces in a natural way the observed "strong alignment of 
the [brightest] galaxy major axis with the long axis of the parent cluster" (Trevese, 
Cirimele, Flin 1992), already mentioned in the Introduction. 

In order to estimate the oscillatory periods of a galaxy around its equilibrium 
position, the numerical values of the TF components in a realistic cluster are 
needed. From eq. (|63|) one obtains 



AT i:j = -27rGp c (0)(l - e 21 )(l - e^ATy. (64) 

The evaluation of the three integrals Ti, T 2 , and T3 involves elliptic functions; in the 
degenerate case of a rotation ellipsoid these functions can be expressed in terms of 
elementary functions. The dimensionless quantities T, are given in Appendix. We 
obtain p c (0) in eq. ( |6~4| ) using the Virial Theorem for the King density distribution 
(King 1972, BT87 p. 228), from which p c (0) ~ 9cr 2 /4vrGaf . The time required for 
a complete oscillation of the galaxy around each of the angular coordinates is then 
given by 



P osc = \ 77— 77a Posc(u, v,e 2 i,e 3 i) ~ 72 ai ' 250 p ogc ( u , v, 621,631) x 10 8 yr, (65) 
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where aiooo is the central (1-dimensional) cluster velocity dispersion normalized 
to 1000 km s _1 (a number of order unity for rich clusters), and 01,250 is the cluster 
core major axis normalized to 250 kpc. Following Valluri (1993) and Girardi et 
al. (1995), we adopt 01,250 = 1. The explicit expressions for the three P osc can be 
obtained from eqs. ((T^)-(jl^). The dimensionality of the parameter space is too 
large to allow a complete exploration, so we limit ourselves to a specific case. For 
example, let us assume 02/01 = 0.8, 03/01 = 0.6, and oiooo = 1- For the galaxy 
model, we also adopt 02/01 = 0.8 and 03 ja\ = 0.6, corresponding to an E2 and E4 
galaxy when the model is projected on the x\ — X2 and x\ — X3 planes, respectively. 
With this choice of parameters, u = I1/I3 = 25/41 and v = I2/I3 = 34/41, and 
the three values of P osc are (2.7,1.6,3.9) x 10 9 yrs, for motions around yj, 1?, ^, 
respectively. So, a comparison with the mean stellar orbital times obtained in the 
previous section (see Fig. 4) shows that in the outer halo of giant Es the stellar 
orbital times can be of the same order of magnitude as the oscillatory periods of 
the galaxies themselves. It could be of some interest a numerical investigation of 
the final fate of these halo stars: in case of escape this is a collisionless evaporation. 
Note that, at variance with the collisional evaporation of globular clusters, in this 
case no mass segregation is expected, since the effect discussed here is independent 
of stellar mass. On the contrary, for the bulk of stars in galaxies the orbital times 
are considerably shorter than P sc> and this is the reason why in the N-body 
simulations of CD94 galaxies are found to behave like rigid bodies. 

4.3 Triaxial Galaxies in Circular Orbit in Spherical Clusters 

As shown in §3, in this case there are two different equilibrium configurations. The 
first equilibrium position, with the galaxy's major axis directed toward the galaxy 
center, reproduces the main observational results already discussed in the Intro- 
duction. In the second equilibrium position the galaxy major axis is perpendicular 
to the orbital plane, as observed for a small number of cases. 

In order to analyze the characteristic frequencies around the equilibrium posi- 
tions in this triaxial galaxy described in §4.1 is assumed on circular orbit 
inside a spherical cluster, described by a Hernquist (1990) density distribution 



[eq. (53) with 7 = 1, a\ = 02 = 03 = r c ]. This assumption about the mass 
profile follows from the most recent results of cosmological high-resolution colli- 
sionless N-body simulations (see, e.g., White 1996, and references therein). From 
the Virial Theorem applied to this density law (Hernquist 1990), one finds that 
18r c <7 2 = GM C , where a is the virial 1-D velocity dispersion. 

The radial trend of q(r) is plotted in Fig. 5 (upper panel, solid line): q mono- 
tonically decreases moving outwards from the cluster center^, where q = 2/3. In 
the same figure, we also plot the normalized angular velocity £1 = Q/^GM c /r% 



For 7- models, q(r) — (1 — y/3)/(l + r/r c ). 
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Fig. 5. Upper panel: the radial trend of q (solid line) and Q (dotted line) for the Hernquist 
model. The radius is normalized to the cluster core radius. Lower panel: the radial trend of the 
three characteristic P oac /Q. The solid line corresponds to the motion around ip, and the dashed 
and dotted lines to P^ c and P<J C respectively. 



where, for each angle, P osc (u,v,q) follows from eqs. (58)-(59), and is plotted in 
Fig. 5 (lower panel) as a function of r/r c for the adopted u and v. Note how P osc 
steadily increases moving outward in the cluster. As a final result there exists a 
radius approximately of the order of the cluster core where the oscillatory times of 
galaxies are comparable with their mean stellar orbital times, and so it is possible 
that some stars will be affected by resonances. The fate of these orbits can be 
investigated by numerical integration of the equations of motion. 



In this paper we have presented a discussion of the motion, near the equilibri- 
um configurations, of a rigid body subject to a tidal field. A complete analysis 
of the equilibrium points when the center of mass of the body is at rest in the 
field generated by a generic mass distribution, and when it is placed on a circu- 
lar orbit inside a spherically symmetric potential, was given. The conditions for 
stable equilibria and the frequencies of oscillations around such positions are ana- 
lytically derived. As an astrophysical application of the previous results, we have 
discussed the observed alignment of galaxies in clusters, concluding the numerical 
investigations presented in CD94. In particular, we have shown that the radial- 
ly aligned configurations found through N-body simulations by CD94 are indeed 
equilibrium positions. We have also shown the existence of an equilibrium configu- 
ration not found by CD94, i.e., the case of a galaxy in circular orbit in an external 



(dashed line). The equation analogous to eq. ( p5| ) i 



is now 




(66) 



5 Conclusions 
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spherical potential, whose major axis is perpendicular to the orbital plane; such 
configurations are observed for a minority of galaxies in real clusters. Moreover, 
comparing the orbital times of stars inside the model galaxies with the oscillatory 
times around the equilibrium configurations of the galaxies in the cluster tidal 
field, we found that for realistic parameters, valid for the majority of the galax- 
ies, the stellar orbital times are much shorter, and so the orbits are adiabatically 
invariant. This could be an explanation of the result numerically found by CD94 
that N-body galaxies behaves in the cluster TF as rigid bodies. 

Finally, we have pointed out the possible effect of resonances between stellar 
orbits in the external parts of the galaxies, and the oscillations frequencies induced 
on the galaxies by the tidal field of the parent cluster. We found that the induced 
frequencies on galaxies can be of the order of the stellar orbital times in their 
outer halos, and this fact suggests the existence of a possible resonant collisionless 
evaporation of stars from galaxies in clusters. This mechanism could explain the 
recently observed population of intracluster stars (see, e.g., Mendez et al. 1997, 
and references therein). 
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Appendix 

We start writing explicitely the expressions for the dimensionless functions Tj, 
defined as Tj = — 2nG p c (ti) {0,20,3 / a\ )Tj [see eq. (|63|)1 in the general triaxial case, 
defined by a\ > 02 > 03. 

From Byrd &; Friedman (1971, hereafter BF71), defining 

e = arccos ( ^) ; k 2 = 1 -4/ a i (67) 

\aiJ 1-ayaf 



one obtains 

~ = 2[F(9,k)-E(9,k)} 
1 k 2 sm 3 (9) 

(BF71, 238.03-310.02), 



(68) 



~ 2[E(9, k) - (1 - k 2 )F(9, k) - (a 3 /Q 2 )fc 2 sin(g)] 
2 A; 2 (l-A; 2 )sin 3 (0) ' 1 ' 
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(BF71, 238.04-318.05), and 



T 3 



2[(a 2 /a 3 ) S m(e)-E(9,k)] 
(1 - k 2 ) sin 3 (#) 



(70) 



(BF71, 238.05-316.02), where F{9, k) and E{9, k) are the Elliptic Integrals of the 
first and second kind, respectively (see, e.g., BF71, 110.02- 100.03). 

The first degenerate case is that of a prolate cluster density, i.e., a\ > a 2 = a 3 . 
In this case, the ellipsoidal eccentricity is given by 



'1- 



(71) 



and the two different components of the TF are given by: 



Ti = 



a\ja\ 



1-e 2 , fl + e 
m 



2e 3 



1-e 



1-e 2 



(72) 



(Gradshteyn & Ryzhik 1965, eq. 2.248.1; hereafter GR65), and 



e 2 2e 3 (l - e 2 ) 



In 



1 + e 
1-e 



(73) 



(GR65, 2.248.4). The other degenerate case is the oblate one. We assume now 
o-i = CL2 > 0-3, and so 



e = yi- 

The components of the TF are readily found: 

1 



,2 4 



Ti = T 2 = 



a 3 /oi 



Vl — e 2 arcsin(e) 1 — e 2 



(GR65 2.248.4), and 



T 3 = 



a 3 /ai 



1 y/1 — e 2 arcsin(e) 



(GR65 2.248.1). 



(74) 



(75) 



(76) 
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